Non-equilibrium thermodynamic study of magnetization dynamics in the presence of 

spin-transfer torque 
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The dynamics of magnetization in the presence of spin-transfer torque was studied. We derived 
the equation for the motion of magnetization in the presence of a spin current by using the local 
equilibrium assumption in non-equilibrium thermodynamics. We show that, in the resultant equa- 
tion, the ratio of the Gilbert damping constant, a, and the coefficient, /3, of the current-induced 
torque, called non-adiabatic torque, depends on the relaxation time of the fluctuating fleld Tc. The 
equality a = [3 holds when Tc is very short compared to the time scale of magnetization dynamics. 
We apply our theory to current-induced magnetization reversal in magnetic multilayers and show 
that the switching time is a decreasing function of Tc. 
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Spin-transfer torque-induced magnetization dynamics 
such as current-induced magnetization reversal P, 0, H j 
domain wall motion [3] , and microwave generation [^| 
have attracted a great deal of attention because of their 
potential applications to future nano-spinelectronic de- 
vices. In the absence of spin-transfer torque, magnetiza- 
tion dynamics is described by either the Landau-Lifshitz 
(LL) equation Q or the Landau-Lifshitz-Gilbert (LLG) 
equation 01 . It is known that the LL and LLG equations 
become equivalent through rescaling of the gyromagnetic 
ratio. 

However, this is not the case in the presence of spin- 
transfer torque. For domain wall dynamics, the following 
LLG-type equation has been studied by several groups 

dt{M) + V ■ V(M) =-fH X (M) 
+ ^(M) X dt(M) + ^(M) X [(,; . V) (M)] , (1) 

where M represents the magnetization, v is the velocity, 
7 is the gyromagnetic ratio and a is the Gilbert damping 
constant. The second term on the left-hand side repre- 
sents the adiabatic contribution of spin-transfer torque. 
The first and the second terms on the right-hand side are 
the torque due to the effective magnetic field H and the 
Gilbert damping. The last term on the right-hand side 
of Eq. ([J) represents the current-induced torque, called 
"non-adiabatic torque" or simply the (3 term. The direc- 
tions of the adiabatic contribution of spin-transfer torque 
and non-adiabatic torque are shown in Fig. [T] (a). 

As shown by Thiaville et ai, the value of the coefficient 
f3 strongly influences the motion of the domain wall 8] . 
However, the value of the coefficient f3 is still controver- 
sial, and different conclusions have been drawn from dif- 
ferent approaches [9, 1^, im, H H H, [H ■ For example, 
Barnes and Maekawa showed that the value of /3 should 
be equal to that of the Gilbert damping constant a to sat- 
isfy the requirement that the relaxation should cease at 
the minimum of electrostatic energy, even under particle 
flow. Kohno et al. performed microscopic calculations 
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FIG. 1: (a) The direction of the magnetization M , the adia- 
batic contribution of spin-transfer torque, {v ■ V) M , and the 
l3 term, M x [{v ■ V) M], are shown. The direction of the ve- 
locity V is indicated by the dotted arrow, (b) The magnetic 
multilayers, in which the pinned and the free layers are sepa- 
rated by a nonmagnetic spacer layer are schematically shown. 
The magnetization vectors of the pinned and free layers are 
represented by Si and Ss, respectively. The effective mag- 
netic field to which S2 is subject is represented by H . (c) 
The direction of the magnetization of the free layer, S2, the 
spin-transfer torque {S2 x Si) x S2, and the non-adiabatic 
torque, S2 x Si, are shown. The direction of Si is indicated 
by the dotted arrow. 



of spin torques in disordered ferromagnets and showed 
that the a and /3 terms arise from the spin relaxation 
processes and that a 7^ /3 in general • Tserkovnyak et 
al. [m derived the /3 term using a quasiparticle approx- 
imation and showed that a = (3 within a self-consistent 
picture based on the local density approximation. 

In the current-induced magnetization dynamics in the 
magnetic multilayers shown in Fig. [1] (b) [3, 17 , 3| j the 



non-adiabatic torque exerts a strong effect, and therefore 
affects the direct-current voltage of the spin torque diode, 
as shown in Refs. 17, [l^l- The magnetization dynam- 
ics of the free layer, S2, has been studied by using the 
following LLG-type equation, 

dtS2 - -gh {S2 X ^i) xS2^jHxS2 + -^82 x dtS2 
e b2 



VIS2 X Si, 



(2) 



where / is the charge current density, g is the anrplitude 
of the spin torque introduced by Slonczewski |1|, ft is 
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the Dirac constant and r] represents the magnitude of 
the "non-adiabatic torque" which is sometimes called the 
field-like torque [13, [Hi- 

In this paper, we study the magnetization dynamics 
induced by spin-transfer torque in the framework of non- 
equilibrium thermodynamics. We derive the equation of 
motion of the magnetization in the presence of a spin 
current by using the local equilibrium assumption. In the 
resultant equation, the Gilbert damping term and the /3 
term are expressed as memory terms with the relaxation 
time of the fluctuating field Tc. We show that the value 
of the coefficient (3 is not equal to that of the Gilbert 
damping constant a in general. However, we also show 
that the equality a = (3 holds if Tc <C 1/ (7^^) ■ We apply 
our theory to the current-induced magnetization reversal 
in magnetic multilayers and show that the switching time 
is a decreasing function of Tc- 

Let us first briefly introduce the non-equilibrium sta- 
tistical theory of magnetization dynamics in the absence 
of spin current [19|]. The LLG equation describing the 
motion of magnetization M under an effective magnetic 
field H is given by 

dtM ^jH X M + xdtM. (3) 

The equivalent LL equation is expressed as 

dtM = ^ ^ ^ ifxM- ,^,"^ Mx{MxH). (4) 
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The Langevin equations leading to Eqs. ([3]) and ([4]) by 
taking the ensemble average of magnetization m, are 



dtm = 7 Htot X rn 



dt5H: 



{5H -Xsm)+Rit), 



(5) 
(6) 



where the total magnetic field J?tot is the sum of the 
effective magnetic field H and the fluctuating magnetic 
field SH and Xs is the susceptibility of the local magnetic 
field induced at the position of the spin. According to Eq. 
([6]) the fluctuating magnetic field SH relaxes toward the 
reaction field XsItt- with the relaxation time Tc- The ran- 
dom field R{t) satisfies {R{t)) — and the fluctuation- 
dissipation relation, {Ri{t)Rj{t')) = ■^XskBT6ij6{t — t'), 
where fee is the Boltzmann constant, T is the tempera- 
ture, (• ■ • ) denotes the ensemble average, and i,j = 1,2,3 
represents the Cartesian components. It was shown that 
Eqs. dl]) and ([6]) lead to Kawabata's extended Landau- 
Lifshitz equation (20| derived by the projection operator 
method [l9|. In the Markovian limit, i.e., Tc <C l/{'-fH), 
we can obtain the LLG equation ([3]) and the correspond- 
ing LL equation ^ with a = jTcXsM [l^. 

In order to consider the flow of spins, i.e., spin cur- 
rent, we introduce the positional dependence. Since we 
are interested in the average motion, it is convenient to 
introduce the mean velocity of the carrier, v. The average 



magnetization, {m{x,t)), is obtained by introducing the 
positional dependence and taking the ensemble average 
of Eq. ([U . In terms of the mean velocity, the ensemble 
average of the left-hand side of Eq. ([5]) leads to 



dt{m) + iv-V){m). 



(7) 



Assuming {5H x m) « {SH) x (m), which is applicable 
when the thermal fluctuation is small compared to the 
mean value, we obtain 

dt{7n) + [vV] {m)=-i{Htot{x,t)) x (m(a;,t)}. (8) 

The mean magnetization density is expressed as 
{M{x,t)) = p{x,t){m{x,t)) , i.e., by the product of the 
scalar and vectorial components both of which depend 
on the position of the spin carrier at time t. The spin 
carrier density satisfies the continuity equation. 



dtp{x,i) + V ■ {vp{x,t)) = Q. 



(9) 



By multiplying the left-hand side of Eq. ^ by p{x,t) 
and using the continuity equation ([9]), the closed expres- 
sion for the mean magnetization is obtained as [2l[ 

p{dt{m) + V ■ V(m)) = dtp{m) + (m)V ■ vp + pv ■ V(m} 
^dt{M) +Divv{M), (10) 

where Divv{M) is defined by 
Divv{M)=J2^^r^= {M){V-v) + {v-V){M).{ll) 

By multiplying the right-hand side of Eq. ([8]) by p{x,t) 
and using Eq. (fTU]) . we obtain 



dt{M) +T)iYv{M) ^-f{H+{6H)) x (M). (12) 

Equation (fT2|) takes the standard form of a time- 
evolution equation for extensive thermodynamical vari- 
ables under flow 2l|. The average of Eq. © with the 
positional dependence is given by 



dt{6Hix,t)) 



[{SH{x,t))-x{Mixit),t))], (13) 



where x{t) is the mean position at time t of the spin car- 
rier, which flows with velocity v = dtx(t) and x — Xs/ P 
is assumed to be a constant independent of the position. 
Equations (fT^ and constitute the basis for the sub- 
sequent study of magnetization dynamics in the presence 
of spin-transfer torque. 

The formal solution of Eq. (|13p is expressed as 



/t 
i;{t-t'){M{x{t'),t'))dt', (14) 
-00 

where the memory kernel is given by i^{t) = cxp[— t/Tc]. 
Using partial integration, we obtain 



{5H{x,t))^x{M) 



ij{t^t')x{M{t'))dt', (15) 
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where the exphcit expression for M{t) = 7Vf(a;(i),t) is 
given by the convective derivative, 

M{t) = dtM{x{t), t) + {v V) M{x{t),t). (16) 

Substituting Eq. P3|) into Eq. we obtain the equa- 

tion of motion for the mean magnetization density, 

dt{M) + Div'i;(M) =-fH x (M) 
+ 7 / dt'ijit - t')x{M{t)) X {M{t')). (17) 

^ — oo 

Equation (|17p supplemented by Eq. (|16p is the principal 
result of this paper. 

When the relaxation time of the fluctuating field, Tc, is 
very short compared to the time scale of the magnetiza- 
tion dynamics, the memory kernel is decoupled and Eq. 
(|17p can be written in the form of an LLG-type equation 
as 

dt{M)+'Divv{M) =jHx{M) + -^{M)x{M), (18) 

where a = ^TcxM is the Gilbert damping constant. Sub- 
stituting the explicit form of the convective derivative, 
Eq. (HI]), into Eq. ^ and using Eq.dTTl) we obtain the 
following LLG-type equation: 

dt{M) + {M){V-v) + {v-V){M) =jHx (M) 

+ ^(M) X d,{M) + ^{M)x [{v ■ V) (M)] . (19) 

If V • i; = 0, Eq. (HH) reduces to Eq. (14) of Ref. 0, 
which is derived by replacing the time derivative of mag- 
netization dtM on both sides of the LLG equation ^ 
by the convective derivative dtM + v ■'V ■ M. The term 
(M}(V-v) appears not on the right-hand side of Eq. (fT9|) 
but on the left-hand side, which means we cannot obtain 
Eq. (|19p using the same procedure used in Ref. 3. As 
shown in Refs. [1 [H, Eq. ^ with (M) (V ■ v) = 
leads to a steady-state solution in the comoving frame, 
{M{t)) = {Mo{x - vt)), where {Mo{x)) denotes the 
stationary solution in the absence of domain wall mo- 
tion. However, if (M) (V-v) 7^ 0, the steady-state so- 
lution may break the Galilean invariance. The situa- 
tion (M) (V ■ v) ^ can be realized, for example, in 
magnetic semiconductors [2^ 24 1 , where the spin carrier 
density is spatially inhomogeneous, i.e., Vp ^ 0. 

The last term of Eq. represents the non-adiabatic 
component of the current-induced torque, which is also 
known as the "/3 term". By comparing Eq. (|19p with 
Eq. ([T|), one can see that the coefficient of the last term 
is equal to the Gilbert damping constant a. However, Eq. 
(fT9)) is valid when the relaxation time of the fluctuating 
field, Tc, is very short compared to the time scale of the 
magnetization dynamics. It should be noted that the 
general form of the equation describing the magnetization 
dynamics is given by Eq. p7p where the last term on the 



right-hand side is the origin of the a and /3 terms. It is 
possible to project the torque represented by the memory 
function onto the direction of the a and (3 terms. This 
projection leads to a 7^ /3 in general. 

In order to observe the effect of Tc on the magneti- 
zation dynamics we applied our theory to the current- 
induced magnetization switching in the magnetic multi- 
layer shown in Fig[T] (b). We assumed that the fixed and 
free layers are single-domain magnetic layers acting as a 
large spin characterized by the total magnetization vec- 
tor defined as 5'^ = / dV{Mi), where i = 1(2) for the 
fixed (free) layer and J dV denotes the volume integra- 
tion over the fixed (free) layer. Both the magnetization 
vector of the fixed layer Si and the effective magnetic 
field, H, acting on the free layer lie in the plane. 

Integrating Eqs. (fT2|) and (fT3|) over the volume of the 
free layer, we obtain the equations, 

dtS2 + J dSh-J^j{H+ {5H)) x S2, (20) 

dt{SH) = --{{dH)-xvS2), (21) 

Tc 

where J — v (g) {M) is the spin current tensor / dS rep- 
resents the surface integration over the free layer, n is 
the unit normal vector of the surface, and xv = x/^ is 
defined by the volume of the free layer V. 

The same procedure used to derive Eq. pT]) yields 

dtS2 + [ dSn ■J = jH X S2 



+ 7 / dt'^it - t')xvS2{t) X dt'S2{t'), (22) 



where ipit) — exp[— ^/tc]. 

When the relaxation time of the fluctuating field is 
short compared to the time scale of magnetization dy- 
namics, the LLG-type equation in the presence of the 
spin-transfer torque is obtained as 

dtS2 + J dSn-J = jH XS2 + ^S2 X dtS2, (23) 

where a = JTcXvS2- By introducing the conventional 
form of the spin-transfer torque [l| , we obtain the follow- 
ing LLG-type equation: 

dtS2--ghiS2xSi)xS2=-/HxS2 + -^S2xdtS2. (24) 



However, Eq. ([M)) is valid only when Tc < l/{-fH). As 
mentioned before, the torque represented by using the 
memory function generally has a component parallel to 
the non-adiabatic torque. In order to observe the ef- 
fect of the non-adiabatic torque induced by the memory 
function on the magnetization dynamics, we performed 
numerical simulation using Eqs. pp]) and ([?!]) . 
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FIG. 2: The z-component of the magnetization •S'2 is plotted 
as a function of time for various values of Tc. The initial 
direction of the free layer is taken to lie in the direction of 
the effective magnetic field, which is aligned to the z axis. 
The initial angle between Si and S'2 is taken to be 45°. The 
Gilbert damping constant a. is 0.01. 

For the simulation, we used the foUowing conditions. 
At the initial time of t = 0, we assumed that the mag- 
netization of the free layer is aligned parallel to the ef- 
fective magnetic field H and the angle between the mag- 
netizations of the fixed and the free layers is 45°. This 
arrangement corresponds to the recent experiment on a 
magnetic tunnel junction system We also assumed 
that the fluctuation field has zero mean value at i = 0, 
i.e., {SH{0)) = 0. 

In Fig. [21 we plot the time dependence of the z com- 
ponent of the magnetization of the free layer, S'2, under 
the large-enough spin current to flip the magnetization 
of the free layer, IghSlSi/{ea^H) = -10. The value of 
Tc is varied while the value of a = 0.01 is maintained. 
The solid, dotted, and dot-dashed lines correspond to 
^Htc — 0.1, 1.0, and 10.0, respectively. As shown in 
Fig. [21 the time required for the magnetization of the 
free layer to flip decreases with increasing Tc, which can 
be understood by considering the non-adiabatic torque 
induced by the spin current. The non-adiabatic torque 
induced by the spin current is obtained by projecting 
the torque given by the last term of Eq. (|22p onto the 
direction of ^2 x Si, which results in the positive con- 
tribution to the spin-flip motion of ^2. Since the last 
term of Eq. (|22p includes a memory function, the non- 
adiabatic torque induced by the spin current increases 
with increasing Tc. Therefore, the time required for S'2 
to flip decreases with increasing Tc. For ^Htc > 10 we 
observe no further decrease of the time required for S2 
to flip because the memory function is an integral of the 
vector S2{t) x dt'S2{t') and the contributions from the 
memory at i — i' ^ l/{'jH) is eliminated. 

In conclusion, we derived the equation for the motion 
of magnetization in the presence of a spin current by us- 
ing the local equilibrium assumption in non-equilibrium 



thermodynamics. We demonstrated that the value of the 
coefficient (3 is not equal to that of the Gilbert damping 
constant a in general. However, we also show that the 
equality a — P holds if Tc l/{'yH). We then apphed 
our theory to current-induced magnetization reversal in 
magnetic multilayers and showed that the switching time 
is a decreasing function of Tc. 
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